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Abstract. Mechanics of fluid membranes may be described in terms of the concepts 
of mechanical deformations and stresses, or in terms of mechanical free-energy 
functions. In this paper, each of the two descriptions is developed by viewing a 
membrane from two perspectives: a microscopic perspective, in which the membrane 
appears as a thin layer of finite thickness and with highly inhomogeneous material 
and force distributions in its transverse direction, and an effective, two-dimensional 
perspective, in which the membrane is treated as an infinitely thin surface, with 
effective material and mechanical properties. A connection between these two 
perspectives is then established. Moreover, the functional dependence of the variation 
in the mechanical free energy of the membrane on its mechanical deformations is first 
studied in the microscopic perspective. The result is then used to examine to what 
extent different, effective mechanical stresses and forces can be derived from a given, 
effective functional of the mechanical free energy. 



PACS numbers: 68.15.+e, 83.10.-y, 87.16.Dg 
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1. Introduction 

Mechanics of fluid, lipid-bilayer based membranes has been one of the most important 
topics in membrane physics during the past three decades Pfl |2j. It contains much 
richer physics than that of conventional fluid-fluid interfaces, due to the fact that the 
concept of surface tensions is, under most physically relevant situations, not sufficient 
to describe mechanical states of the membranes, unlike in the case of conventional fluid- 
fluid interfaces. Additional mechanical quantities to describe their resistance towards 
bending deformation have to be included in the description as well J3H]. Understanding 
of the equilibrium shapes of red-blood cells [Sj and giant, artificial lipid-bilayer vesicles 
[HI EH HI IE] relies on a well-defined description of the mechanical properties of the relevant 
membranes, so do the applications of membrane mechanics in experimental techniques 
such as the micro-mechanical manipulation technique (HI E]j and the sensitive force- 
measurement technique based on the mechanics of cells or vesicles [TT1 H2 EH] . Recently, 
descriptions of membrane mechanics have also been employed in the investigation of 
surface-induced forces between membrane-adsorbed or membrane-embedded colloidal 
particles [HI HE]. 

Descriptions of mechanics of fluid membranes are most often formulated from 
the point of view of elastic shell theory, or an effective, two dimensional perspective, 
where the membrane is treated as an infinitely thin surface. A description is given 
either in terms of surface stresses and bending moments, expressed as functions of 
the deformations, or in terms of mechanical work or free energy associated with 
the deformations Ei]- A state of mechanical equilibrium of the membrane, or 
its equilibrium shape, can then be determined by setting up equations of balance 
for the stresses and for the bending moments, if they are given. Obviously those 
equations correspond to the conservations of linear momentum and angular momentum, 
respectively PJ llfij . Alternatively, the problem of finding equilibrium shapes of the 
membrane can be approached, more directly, by applying a variational principle to derive 
the equations of mechanical equilibrium from the mechanical free energy [T? [ l6| 1X8 } 18]. 
Previous works have shown that these two approaches can be reconciled with each other 
|19| l2*U]. In particular, the work presented in [20] based on the use of Nother's theorem 
provides a very general, albeit somewhat abstract, argument regarding why the two 
approaches are equivalent. 

One of the goals of this paper is to give a more mechanical - therefore, more 
experimentally relevant - account than that given in [20], of why the two approaches 
of obtaining equations of membrane mechanical equilibrium are equivalent. For this 
purpose, general descriptions of the mechanics and thermodynamics of a fluid membrane 
will be developed by viewing the membrane from two different perspectives. One is a 
microscopically realistic perspective, depicted in Figure H^a), where the membrane has 
a highly inhomogeneous distribution of matter in its transverse direction. This will 
often be referred to as the microscopically viewed system or the microscopic model. 
The other perspective is an effective, or idealized, one, illustrated in Figure d[b), where 
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Figure 1. A schematic sketch of the membrane- fluid system under our consideration. 
Part (a) depicts the interfacial region in the real system, showing a membrane 
composed of bilayer-forming lipids and transmembrane proteins in contact with two 
fluids containing small solutes. Part (b) illustrates the representation of the interfacial 
region in the corresponding Gibbs model system. 



the membrane is modeled as an infinitely thin surface, which will be called the dividing 
surface, with the bulk solvent filling all of the space on the two sides of the surface. 
This model will also be called the Gibbs system or Gibbs picture, as effective, excess 
mechanical quantities will be assigned to this surface following the idea of Gibbs' |21j . 
Those excess quantities are defined by the constraints of mechanical equivalence between 
the two different perspectives. In other words, the total mechanical forces and torques 
obtained from integrating the linear stress and the bending moment over a surface which 
traverses the membrane surface and has a sufficiently large transverse dimension must be 
the same in the two perspectives. If the transverse extensions of the surface of integration 
into the two solvent regions on the two sides of the membrane are denoted e + and e~, 
the term "sufficiently large transverse dimension" refers to the following requirement: 
that e + and e~ be judiciously chosen such that at these distances from the dividing 
surface the thermomechanical properties of the solvent must be indistinguishable from 
those of a homogeneous bulk solvent under the same thermodynamic conditions. 

A second goal of the paper is to add to the canonical descriptions [TJ HB1 ED] 
descriptions of physical situations where applied external forces that act, or do work, on 
membranes are not localized on the corresponding dividing surfaces. An experimentally 
relevant example of such situations is the following |2*2*1 |23] : an external electric field is 
applied to a membrane system, induces a dipole moment in the membrane, and through 
its coupling to the dipole moment exerts a mechanical torque on the membrane. The 
work of such a torque will depend on membrane deformations and therefore, contribute 
to the mechanics of the membrane. Such situations have not been considered in the 
earlier works ^ El 120] • 

A third goal of the paper is to address an issue that arises from the exercise of 
deriving mechanical stresses from a given thermomechanical free-energy functional. It 
turns out that mechanical stresses obtained by this way inevitably contain in their 
expressions some degree of arbitrariness. A similar problem has a long history in classical 
field theory [24J. In this paper, the origin of the arbitrariness is investigated and the issue 
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of to what extent this arbitrariness can be eliminated in physically meaningful or relevant 
ways is discussed. This part of our study should facilitate the task of determining states 
of mechanical equilibrium of membranes. 

It should be noted that even though we use the word membrane throughout the 
paper, then the theory is general enough also to apply to for instance the interface of 
two coexisting fluids. However, since a membrane is not necessarily situated at such an 
interface, it can for instance be freely floating, and since the prime example that we have 
in mind is a lipid-bilayer based membrane, we will use the word membrane throughout 
the text. 

The outline of the paper will be as follows. In Section El the Gibbs description of 
the mechanics of fluid membranes will be introduced. In Section |3] a description based 
on the corresponding microscopic view of the membrane systems is established and is 
then connected to the Gibbs description. This connection is then used in Section |U 
to elucidate how the mechanical stresses in the Gibbs description are related to the 
functional variation of a general mechanical free energy under shape deformations and 
to identify the issue that the mechanical stresses derived from the free-energy variation 
are not uniquely determined. This issue is then addressed in Section A conclusion 
is given in Section and it is finally followed by two appendices. The first of these, 
Appendix |Appendix A[ is a brief summary of the differential geometry of surfaces. It is 
included mostly to have a place to consult for the notation and conventions used. 

2. The Gibbs description of membrane mechanics 

In this section we will develop a description of membrane mechanics based on viewing 
from an effective, or idealized, perspective a system of a membrane together with the 
fluids that surround it. In this perspective, the whole system is approximated by a 
dividing surface of zero thickness together with two regions of bulk fluids, the behavior 
of which are assumed not to be affected by the presence of the membrane. We refer to 
a membrane system modeled in this way as the Gibbs model or the Gibbs system. The 
mechanics of the Gibbs system is then described by a set of effective surface stresses - 
defined on the dividing surface only - together with those known of the bulk fluids. We 
will first focus on developing the concept of effective surface stresses, and then give a 
full description of the mechanics of the Gibbs system. In correspondence with laws of 
conservation of linear and angular momenta fundamental to mechanics, we will classify 
mechanical stresses into linear stress tensors and angular stress tensors. We will develop 
the concept of the linear stress tensors from the more fundamental concept of forces, 
and similarly, the concept of the angular stress tensors from the concept of torques. 

2.1. Forces and linear stresses 

Similar to mechanical stress tensors in conventional material systems, linear stresses in 
a membrane are employed to describe the forces acting on an element of the membrane 



Descriptions of membrane mechanics. 



5 



by the rest of the membrane with which it is in contact. The nature of such contact 
forces is that they are boundary forces. Denoting the element of the membrane under 
consideration E, we can define a vector, TV) = T( l /)(^ 1 , £ 2 ), to be the force per unit 
length acting on the boundary of E. The vector i> = v a t a represents the unit vector 
which points in the outward normal direction to the boundary element local to point 
(£\£ 2 ) an d which lies tangentially in the membrane plane. Clearly, the total force on 
E arising from the boundary forces, F-£, is then given by integrating TV) along the 
boundary of E, 

F s = f als TV) , (1) 

JdT, 

where ds is the arc length along the boundary <9E. 

Obviously, TV) depends on the orientation of v. This dependence can be shown 
to be 

TV) = T a u a , (2) 

where the two constituting space force vectors, T Q , which also form surface vectors, are 
independent of v. We have for completeness given a derivation of Equation (j2J) using a 
standard argument in |Appendix B| Readers willing to accept the result may skip 
that appendix. 

It should be clear that T a describes the non-trivial local distribution of forces, and 
it will therefore be identified as the linear surface stress (tensor) in the membrane. It 
will often be used later in the form of its components in a decomposition of the following 
form: 

T Q = T a % + T n a n , (3) 

where its tangential components, T a/3 , constitute a surface tensor of rank 2, and its 
transverse components, T" , form a surface vector. 

Using Equation (J2J) we can now rewrite force F^,: 



ds T a v a = / dA D a T a , (4) 
'<9s Js 

where Gauss's law has been used for the second equality sign. From this we see that 
a surface density of force can be associated with the stress tensor through its covariant 
divergence, 

/ = D a T a . (5) 

If / is non-zero in equilibrium, then an "external" force acting on the membrane to 
balance it is necessarily involved. Such an external force can, for example, result from 
the stresses exerted on the membrane by the surrounding bulk fluids. We will return to 
this point in Section El 

With the help of the components of the linear stress tensor, we can decompose, TV), 
the force per unit length acting on a membrane element through one of its boundary- 
curve elements with outward normal u, into three components, 

TV) = 7> + T s (n x w) + Q T n . (6) 
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These three components have direct physical interpretations: T M = T a pv a v^ can be 
identified as the surface tension acting on the boundary element, T s = T a pv a v 1 e l13 can 
be identified as the shear tangential to the boundary, and Qt = T£v a is transverse 
shear normal to the surface. These concepts have already been introduced and used in 
previous seminal work on membrane mechanics [1. 

2.2. Torques and angular stresses 

Following a line of reasoning similar to that sketched in the preceding subsection, we can 
also develop the concept of angular stresses from the concept of torques. Considering 
the same local element, E, of the membrane, with its boundary curve <9£, we can 
define two densities of torque: fl(u), the torque per unit length acting by the part of 
the membrane neighbouring to £ on an element of the boundary curve with v as its 
outward pointing normal vector, and r, the total torque per unit area resulted from 
the torques distributed along the boundary. In the rest of the paper, it will be assumed 
implicitly that any torque referred to is calculated with respect to the origin of the 
global coordinate system unless specified otherwise. We have, therefore, 



dAr = ds fi (l/) . (7) 
We can also establish the following two equations, 

n (I/) = n a u Q , (8) 

and 

r = D a Q a , (9) 

based on the same arguments that have led to Equation (J2J) and Equation (JHJ). Obviously 
the space-surface vector quantity, fl a , is the counterpart of T a in the context of angular 
momentum, and will therefore be called the angular stress tensor. 

To make clear the physical content of the angular stress tensor, we decompose it 
into two parts. One is the angular stress generated by the action of the linear stress, 
which is given by R x T a . The other, defined by 

N a = n a - RxT a , (10) 

represents a contribution which is non-trivial. Its origin lies in the fact that the real 
physical system of a fluid membrane is not an infinitely thin surface and that the 
force distribution along its transverse dimension is inhomogeneous. This point will 
be discussed more in the next Section. We will call N a the internal angular stress 
tensor in reference to its origin. In correspondence with the decomposition of fl a we 
divide the surface density of the torque, r, into two parts, 

r = Rxf + r int . (11) 

The first part describes the contribution from the surface density of the force, and the 
other part, which can easily be rewritten as 

r int =r-Rxf= D a N a + t a xT a } (12) 

will be referred to as the internal torque. 
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Figure 2. A sketch of the cell, or volume element, E. B 1 and B 2 refer to side surfaces 
of the cell. 



2.3. Mechanical stresses in the Gibbs system 

Having developed in a systematic way the concepts of surface linear stress tensor and 
angular stress tensor, we can formulate precise expressions of the corresponding three- 
dimensional stresses that describe the force and torque distributions across the whole 
Gibbs system consisting of the dividing surface and the bulk fluids. 

To do that, we consider a three-dimensional cell E, as shown in Figure 121 which 
includes the local element of the dividing surface E, but which has sufficiently large 
extension in the transverse direction. The four side faces of cell E, which are transverse 
to E are labeled as B± and B±, are parametrized mathematically by 

B« = {R(?,?) + hn(?,?) 

£3-a _ A£ 3 - a /2 < < £ 3 ~ a + A£ 3 ~°72 , 

C = C ±AC/2 , -e-<h<e + }. (13) 

The intersections of B± with the dividing surface, B±, hence constitute the 
boundary, <9E, of of the dividing surface element E. 

To formulate mathematically concise expressions of the relevant mechanics 
quantities of the whole Gibbs system, we introduce two Heaviside step functions, which 
are defined as 6 ± (r) = 8(±<j)(r)). The scalar function, 0(r), appearing in the definition 
is defined to be such that it is zero on the dividing surface, positive in the bulk region 
that n points into and negative in the other. We also use two three-dimensional stress 
tensors, T , to represent the thermomechanical properties of the bulk fluids, which are 
assumed known in terms of their dependences on thermomechanical control variables. 
One of the mechanical quantities, which will be needed soon, is the total force resulted 
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from integrating the three-dimensional stress tensor over any of the four side surfaces, 

F GihhB (B±) = [ ds T% + [ dAv- (6 + T + + rt") , (14) 
Jb<£ Jb^ 

where v is the outward-pointing unit vector normal to the relevant side surface. 
3. A microscopic description of membrane mechanics 

Having introduced different mechanical quantities that constitute the Gibbs description 
of mechanics of a fluid membrane, we will in this Section approach the mechanics of 
the membrane from the microscopic perspective briefly described in the Introduction. 
In this perspective, distributions of forces and torques, or mechanical stress tensors, are 
considered three dimensional, and in particular, the inhomogeneity in the distributions 
along the dimension transverse to the membrane surface is described explicitly. Our 
main goal is to establish a concrete connection between this microscopic description 
and the Gibbs description. The physical concept that underlies this connection is that 
of mechanical equivalence. Specifically speaking, we will build up this connection by 
considering the linear stresses, the corresponding forces, the angular stresses, and the 
corresponding torques one by one. A similar approach, dealing with the cases of the 
stresses, can be found for instance in Ref. [T5] . 



3.1. The linear stress tensors 

In the microscopic description, a linear stress tensor, T real , is employed to describe the 
local distribution of force around any point in the three dimensional space occupied by 
the membrane and its surrounding fluids. It is a tensor of rank 2. To relate this tensor 
to the effective linear stress tensor, T a , which is defined only on the two dimensional 
dividing surface, we apply the concept of mechanical equivalence explicitly. It states that 
in the Gibbs and the microscopic descriptions, the total forces, resulted from integrating 
their respective linear stress tensors over each of any two independent side faces chosen 
out of B± must be the same. In expression, the statement reads as 

/ dA v ■ f real = I ds T% + / dAv- (0+T+ + A"? - ) , (15) 

Jb%_ JBg Jb<£ 

where Equation dHj) has been used. 

Alternatively, we can define a stress tensor, T excess , by 

Texcess = T rea l — + T + — 9 T~ . (16) 

This tensor represents all the excess stress that results from the existence of the 
membrane interface. Using it we can rewrite Equation (fTB^) as 

/ dsT p vp= I dAi?-J excess . (17) 

J B1 JBg 
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To derive from Equation (J 17)) an explicit expression of T a in terms of the 
microscopic stress tensor we must rewrite the two sides of the Equation ([17)1 . For 
the rewriting of the right-hand side we choose to work with the curvilinear coordinates 
(^ 1 5^ 2 ) h), which parametrize the spatial position of any point R in the system by 

R(e,e,h) = R(e,e)+hn(e,e)- m 

We then parametrize any of the boundary- curve elements, B±, by R(l; l (\), £ 2 (A)), in a 
chosen orientation such that the outward pointing normal vector to the corresponding 
side face, B±, is given by 

tv4- {t p (19) 



v 



where 

t p = dpR = {5} - hK\) t 7 . (20) 
It follows then that 

d^ 

dAv = dXdh (tp xn)-^-. (21) 

d\ 

Using £ 7 x n = t a e ai we get 

I dAu- T excess = / dXdh t a ■ T excess £ a7 (81 - hK 1 'A — . (22) 

J B% J B%. 

Rewriting the left-hand side of Equation (JT7)l by use of the parametrization we get 



[ dsT a u a = I d\T a [t a -(tpxn)}§- 



dX T a e a ^ . (23) 



Comparing Equation (|23j) with Equation (j22j) . we finally arrive at the expression that 
we need: 

T a (e,e) = jdh - hL<#{t\?)] t(,(t\?) 

•T excess (£\£ 2 ,/i) . (24) 



3.2. The resultant force densities 

f, the surface density of the force resulted from the surface excess linear stress, has a 
counterpart / excess = V ■ T excess in the microscopic description, which is the volume 
density of the force resulted from the three-dimensional excess linear stress tensor. 
Following the derivation of Equation ()24j) we can now find in a straightforward way 
an expression of / in terms of / exccss . 

Using the definition of /, given in Equation (JHJ), together with Equation (|17|) and 
applying Gauss' law we have 

[ dAf= [ ds T a u a = [ dAu- t excess = [ dV / excess . (25) 
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Working further based on the parametrization defined by Equation (|18|). we can express 
the volume element in the above integration as 

dV = (l - 2hH + h 2 K) dAdh . (26) 

Since Equation (|25|) holds for any E of arbitrarily small area we conclude then that 

/ = Jdh(l-2hH+ h 2 K) f CKCCSS . (27) 

To illustrate the physical meaning of / excess in another way, we derive for it an 
alternative expression. We introduce first the volume density of the resultant force in 
the microscopic description of the system, which is given by 

/real = V ' T real , (28) 

and those in the bulk fluids of the Gibbs model, which are given by 

f = V • T± . (29) 
Using these quantities we can then express / exccss at any point r in the system as follows 

/excess = /real " & + f + ~ ^ f ~ [ dA U ■ (t+ - T~) 6 ( P - R) , (30) 

Jm 

where the subscript M indicates that the integral is over all of the membrane and R is 

the membrane shape function. 

Based on Equation (|5Dj). it is clear that / excess consists of not only the difference 

— — j- 

in the volume force densities, / real and / , but also any non-zero stress difference 
across the dividing surface, n ■ (T + — T~), associated with the Gibbs bulk fluids. For 
example, in the simple case of a conventional system of coexisting fluids separated by 
an interface, where there is no applied external force, the effective linear stress tensor 
T a is characterized by a single mechanical quantity, a homogeneous surface tension a; 

T a = at a (31) 

The effective resultant force / is then 2Ho~n. Equation (|3U|) thus collects all the 

— _ j. 

microscopic effects that contribute to /. If / real = and / =0 are inserted into 

Equation (j3*Uj) as the conditions of mechanical equilibrium for the whole system in the 

microscopic model and for the bulk fluids in the Gibbs model, respectively, Equation (jSUj) 

becomes the familiar statement of mechanical equilibrium for the dividing surface in the 

Gibbs model, namely / must balance the stress (pressure) difference across the dividing 

surface. Where there exists a volume distribution of an external force in the system 

which acts on membrane molecules but not the bulk fluids, such as one that could be 

_ _j_ _ 

induced by an electric field, / real is non-zero and / =0. f excess thus also includes the 
effect arising from any interface- related inhomogeneity in / real . 
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3. 3. Micromechanical expression for the angular stress 

Similar to the linear stress tensor that we have already discussed, the angular stress 
in the realistic microscopic description of a membrane system is a second rank tensor, 
denoted by fi re ai m the following. Although what we develop in this section can be 
applied to cases where the systems under considerations have internal angular momenta, 
we will limit ourselves to considering those where internal angular momenta are not 
relevant. For such a system, the angular stress tensor is related to the linear stress 
tensor in a straightforward way: 

^real = — T real x R , (32) 

where T rea i is assumed symmetric, since a symmetric linear stress tensor can always be 
constructed J. The tensor product is defined as 

(Treal x -R)^- = ^jkl^real,ikRl , (33) 

where the Latin indices range from one to three and indicate the components in 
a Cartesian basis {e^}, and where is a third rank tensor, 6123 = 1 and e^i 
is antisymmetric under any pairwise interchange of indices. Einstein's summation 
convention is also implied here for repeated indices. 

To establish the microscopic origin of the effective excess angular stress defined 
on the dividing surface in the Gibbs description, fi, we use ft rea i as an intermediate 
and follow steps that are very similar to those established in the preceding discussion 
concerning the linear stress tensors. The angular stress tensors in the two Gibbs bulk 
fluids, Q , are related to their linear stress counterparts according to = T ± x R. 
Subtracting these contributions from Q rca i, we obtain a tensor of excess angular stress, 

^excess T excess X R . (34) 

Going a few steps further, we arrive at the following expression: 



/ ds Sfvp = dA Rx (v ■ T excess ) 

JB1 JBZ 



(35) 



Using the parametrizations given in Equation (fTSj) to Equation (|21|). we can re-express 
the left and right-hand sides of the above equations, respectively, as 

ds Q?vp = [ d\ ft% 7 ^V > ( 36 ) 



B1 JB 



d\ 



and 



/ dA R x (p ■ T cxccss ) = / dXdh (R+ hn) x (t a ■ T c 

JB1 JBZ 



■ e ai {5} - hK\) A_ . (37) 

\ In fact, Equation (|32|l . applies even in the case where a system possesses internal angular momentum, 
since it is still possible to construct a symmetric linear stress tensor which takes into account the internal 
angular momentum |26j . 
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Comparing these two expressions we arrive at the final expression which reveals the 
microscopic origin of f2, 

n a = Jdh (g a ? - hL afi ) [R x (tp ■ T cxccss ) + hn x (tp ■ T excess ) ] . (38) 

The two terms in the above expression correspond naturally to two contributions: 
a contribution related to the excess linear stress in the Gibbs model, 

Rx T a = J dh (g a P - hL a ?) Rx (tp - T excess ) , (39) 

and a contribution giving rise to the internal excess angular stress, 

N a = J dhh (g af3 - hL a ?) nx(tp- T excess ) 

= - J dhh (g a ? - hL a ?) t^ 5 (tp ■ t excess • t 5 ) . (40) 

From Equation ([40)1 . it is clear that the primary contribution to N a (the one that is 
least suppressed by the microscopic thickness of the membrane) is the first moment of 
the distribution of the linear stress tensor in the transverse direction. Following the 
nomenclature used in the existing literature on membrane mechanics pQ, we will also 
call N a the bending moment of the membrane. 

Equation (J4*U|) also implies that N a -n = 0. This is a consequence of the assumption 
we have made in our description of the membrane mechanics: that there is no bending 
moment pointing in the normal direction of the membrane surface. This assumption is 
explicitly expressed in Equation ([!%]). and is consistent with our empirical understanding 
of the fluid characteristics of the membrane in its lateral dimensions. This property of 
N a ■ n = allows us to define a related vector quantity, M a , which is given by 

N a = n x M a and M a ■ n = . (41) 
N a and M a have only tangential components, 

M a = M a Hp , N a = N a Hp . (42) 
Expressed in terms of these components Equation ()41)1 assumes two alternative forms, 

M alS = -N^e^ , N a(B = M Q7 e/ . (43) 
Thus we are led to the following micromechanical expression for M a/3 

M al3 = j 'dhh (g a i - hL a ^) (t 7 • T excess • t^) . (44) 

The simplicity of this expression compared to (|4f)|) may already demonstrate why M Q/3 
is often more convenient to use than N a . It is also physically more intuitive, readers 
familiar with Ref. pQ will find that M a/3 are exactly what are called moment resultants 
there. In what follows we will use N a and M a @ interchangeably, choosing the one that 
is the more convenient in the case at hand. 
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3.4- Micromechanical expression for the torque 

A micromechanical expression for the effective area density of excess torque defined in 
Equation Q, t, can also be developed. In the microscopic description we can define 
a volume density of excess torque as the divergence of the excess angular stress tensor, 
i.e., 

''"excess ^ ' ^excess • (45) 

Carrying out steps of derivation similar to those leading to Equation (J27j) for the effective 
area density of excess force, we find that r is related to its counterpart in the microscopic 
description as follows: 

- = / «Mi-™ _ W 

We can further relate r to the force distribution / excess , since 

''"excess R X / excess , (47) 

a consequence of the symmetry of the linear stress tensor in the microscopic description. 
Inserting (|47jl into (j4fij) and performing the integration yields a more revealing expression 
for r, or more specifically, for the effective internal torque, r int , defined in Equation (|T2|). 

T = Rxf + T int , (48) 

where 

Tint = Jdh (1- 2hH + h 2 K) {R-R)x / excess 

= e a %{D 1 M\-n-T a ) . (49) 

It is clear from the above expression that only the tangential components of the 
force distribution, / excess , that are not localized on the dividing surface contributes to the 
effective internal torque. Under the condition of mechanical equilibrium, the tangential 
components of / excess are zero, unless there acts a spatially distributed external force in 
the directions tangential to the dividing surface. Consequently, the internal torque is 
zero under the same condition. 

A couple of remarks on T- mt may be made here. First, Tint has not been part of the 
effective description of membrane mechanics developed in the earlier works [THl IT]. 
There are, however, experimental situations where Ti nt is essential. An example is when 
an electric field induces and couples to an electric dipole moment of the membrane 
[221 E3 Thus we will not neglect T int here. 

Secondly, the fact that the component of T- mt in the direction normal to the dividing 
surface is always zero makes another point clear. Using this fact on Equation (jl2j) by 
setting its normal component to zero leads to the following result: 

T^e aP + M a ^ p K p a = . (50) 

This result states that T a/3 is symmetric only if the bending-moment tensor M a /3 and 
the curvature tensor, K a g, commute, or in other words, if both tensors become diagonal 
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simultaneously, in a properly chosen basis of tangential vectors at every point on the 
surface. This condition may be satisfied in practice in most cases of membrane systems 
but not in general, although its generality has been assumed assumed in some of the 
earlier literature pQ. 



4. Free-energy based description of membrane mechanics 

The descriptions of the mechanics of membranes that we have developed so far are 
formulated directly in terms of the linear and angular stresses, or force and torque 
densities as the divergences of the corresponding stresses. However, often in practice, an 
alternative formulation is used, given in terms of mechanical free energies as functionals 
of deformation of the membranes under consideration. In such cases, the mechanical 
stresses are treated as quantities derived from the free-energy functionals. In this 
Section, we will develop, based on symmetry principles and mathematical argument 
only, a general functional form of the variation of a membrane mechanical free energy 
associated with an arbitrary variation in the deformation of the membrane. We will 
view the membrane mechanics from the microscopic perspective, and classify membrane 
deformation into two categories: one which is of translational nature, and the other 
which corresponds to rotations. From the general functional form of the free-energy 
variation we will then identify those variational quantities that are related to the 
mechanical stresses and indeed, clarify the correct relations between them. At the end 
of the Section, we will apply the general formulation to a more specific, and frequently 
used, example of a mechanical free energy and derive the mechanical quantities such as 



4-1. The mechanical free energy and deformation of a membrane 

In our consideration, we assume that deformations of the membrane occur at constant 
temperatures. Correspondingly, we work with a thermodynamic ensemble where the 
most relevant free-energy density is given by / = e — Ts, where e is the energy density, T 
the temperature and s the entropy density. Following the already established procedure 
for defining the excess of various extensive quantities we define a volume density of the 
excess of the mechanical free energy by 

7exce SS (r-) = / rca i(r) - 6+ (r) /+ (r) - 6~ (r) f- (r) , (51) 

where / rea i(^) is the volume density of the actual free energy contained in the system 
and f ± (r) are defined to be the free-energy densities associated with the "filler" bulk 
fluids in the Gibbs model. 

As before, we consider cell S, whose content of the excess free energy is given by 

-^S, excess = /_ dV /excess • (52) 
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We then define the scalar quantity / to be the excess free energy per unit area associated 
with the dividing surface, i.e. 



We consider / to be the core element in a free-energy based, phenomenological 
model description of the mechanics of a membrane system of interest and take the 
functional dependence of / on the relevant thermomechanical variables as our starting 
point. In the discussion that follows, we will limit ourselves to situations where / is a 
local function of surface density fields of excess molecular numbers n^, where A ranges 
over the different species in the membrane, as well as of derivatives of the shape field, 
R, of the dividing surface. But, we will only consider the functional dependence at a 
generic level, as what is allowed by symmetry as well as mathematical requirements. In 
other words, we will leave the specific functional behaviour of / to be a degree of freedom 
to be fixed in a particular phenomenological model, and only derive consequences from 
the generic functional dependence of / on the variables. 

A few remarks may be given here to clarify our limiting assumption. First, the 
assumption that only derivatives of R, not R itself, are included in / is a direct 
manifestation of our requirement that / be invariant under rigid translations of the 
membrane system. Moreover, we also assume, or require, that / be rotationally 
invariant. In the case where external forces such as gravity or an external electric field do 
work on the system during its deformation, this requirement of the invariances in turn 
means that / represents only the "internal" contributions, or in other words, that the 
contributions associated with the external forces must be considered in addition. This 
requirement is really a matter of choice, made to facilitate the following formulation in 
a technical sense, as it will become clear later on. 

Secondly, additional thermodynamic variables or fields can be included in /, of 
course. As long as these fields remain invariant with respect to rigid translations and 
rotations of the whole system, their addition does not lead to any qualitative changes 
in the arguments that will be developed in the following. As an example, we may 
imagine a situation where the collective orientation of the constituent lipid molecules 
of the membrane under consideration becomes relevant [22j. A vector field, </>, may 
then be used to represent the orientational field. A number of scalar fields which have 
the invariances can then be defined from <p: its length, its projection onto the normal 
direction to the dividing surface and the angle between its tangent component, (4>-t a ) t a 
and the ^-axis in the internal coordinate basis of the local tangent space. 

Lastly, the molecular number density fields ua appearing in / should be understood 
as the excess fields defined as 




(53) 



Immediately we have the following equivalent expression of /, 




(54) 




(55) 
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In other words, the material content of each molecular species in cell E in the Gibbs 
description must be the same as that in the corresponding microscopic description. 

For the purpose of deriving relevant mechanical quantities from the free-energy 
functional /, the variation of / induced by an arbitrary, infinitesimal deformation SR 
in the shape of the dividing surface becomes an essential quantity. Although we are not 
able to evaluate the specific variation without the knowledge of the specific functional 
behaviour of /, we can already make a statement concerning the general functional 
structure of the variation Sf from a purely mathematical point of view: Sf associated 
with the deformation can always be organized into the following form, 

^H^f) = ^H^f)\ T , {V - gnA} = -frs-SR+D a [s a -5R) . (56) 

In the above equation, the vector quantity / rs is regular in the sense that it does not 
involve any differential operators, while 

oo 
n=0 

may involve differential operators. The fixed ^/griA constraint for the variation 
corresponds to the requirement that the excess number of molecules contained in the 
surface element E be conserved for each species under the deformation. 

Two points may already be noted here, concerning Equation (|56|). The first is that 
the variation of the free energy is not sufficient to define S unambiguously. We will 
return to this point later, in Section The second is that / rg = D a S? y This follows 
from the requirement that the free energy be translationally invariant, or equivalently, 
that the right-hand side of Equation (|56|) vanish for constant SR. 

Quantities / rs and S?^" s may be related to the mechanical quantities already 
introduced and discussed, such as /, r, T a . It is the purpose of the rest of this Section 
to develop the connections. In order to do that, we will start from a microscopic 
point of view and consider the three dimensional cell E and an arbitrary, purely 
mechanical deformation SR which it undergoes. Given SR, we can establish, based 
on a mathematical point of view, the following general functional expression for the 
corresponding variation of the free energy, -Fg,excess : 

^£,exccss = - fdV / cxcess -5R+ I dAv-T-SR 

JT, JdT, 

- / dA V ■ [SR - 5(R+ hn)] f excess . (58) 
T in the second integral represents a tensor operator with a general form as 

oo 

= T(n),ijfci...fc n Vfc 1 . . . V fcn , (59) 

n=0 

where quantities ^(n),ijki ...k„ are tensors of different ranks. 

Seen from a physical point of view, the meanings of the three terms in Equation (|58J) 
are almost obvious. The first term represents the mechanical work done by any non-zero 
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volume-distributed external force, which balances / exccss - The second term describes the 
work associated with the forces acting on and distributed over the boundary surface of 
cell E by the environment it is in contact with. The appearance of the tensor operator 
T may, however, have obscured this physical interpretation and we will comment on this 
point a bit later. The presence of the last term has to do with the following fact: that 
the geometry of cell E is defined in such a way that its side boundary surfaces always 
remain orthogonal to the dividing surface. Consequently, there is a free-energy change 
associated with the movement of the matter across the cell boundary induced by SR. 

To explain why we have identified the vector quantity in the first term as / excess , 
let's consider a situation where cell deformation and its derivatives vanish at the cell 
boundary. The variation in the cell free energy must then equal the work done by any 
non-zero external force / rea i,ext under the volume deformation, i.e., 

dV / real , ext SR = - I [dV / real • SR . (60) 
as ^2 

In a similar fashion, the variations in the free energies, f ± , stored in the "filler" bulk 

fluids in the Gibbs description may be expressed as 



» 1 l_dVf Teal 



dV f 

2± 



dV T± : VSR 

52 ^2± 



f dV f° ■ SR± [ dAn-T* ■ SR\^ , (61) 
J2± is 



where the colon means that the two tensors are contracted into a scalar and E ± 
represents the parts of E which are, respectively, above and below the dividing surface 
E. Note that Equation (jfilj) implies that the free energies of the bulk fluids do not 
contain terms that vary with the second spatial derivative of the cell deformation, in 
contrast to the excess free energy assigned to the dividing surface. 

Subtracting Equation (JBTJl from Equation (JoTij) and using Equation (J30J), we find 

^E.excess Le = ~ dV f excess • SR . (62) 



s 



This is exactly the contribution described by the first integral in Equation (JHHJ)- 

Concerning the second integral in Equation (|58|). the connection between its 
physical meaning and its functional structure, or more specifically speaking, the presence 
of the tensor operator, T, may not appear obvious. A reason for the presence of T arises 
from the fact that any given free-energy functional alone does not lead to a unique 
identification of the corresponding stress tensor §. In other words, a kind of "gauge 
freedom" exists for the free-energy based derivation of stress tensors. The expression of 
this gauge freedom necessarily introduces the operator terms contained in Equation (|59|) . 
We will make this point clear by arguing that the non-operator term in T, To, may be 
identified with T excess introduced earlier, which is assumed to have a well-defined physical 
interpretation. The argument may be constructed as follows. 

§ Another reason is that non-zero T(n),ijki...k n f° r n > become necessary if the orientational degrees 
of freedom of the membrane under consideration are relevant, as in the systems of liquid crystals |28j 
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We first consider a rigid translation of the whole cell S, i.e., SR = 5R = C, where 
C is a constant vector. It follows immediately that 

= SF^ exccss = -C J_^dV (/ excess - V • f (0) ) • (63) 

Since S is arbitrary and can also be made arbitrarily small within the limit of a 
continuum description, this means that 

/excess - V • t (0) = V • (T cxcess - T (0) ) = . (64) 

We can then conclude based on a standard theorem in calculus that the difference 
between T exce ss and T( ) is the curl of another tensor V, i.e. 

(65) 

However, as we will see shortly, this difference can be eliminated by a gauge 
transformation to a new tensor T^. In other words, given the knowledge of T excesSi jj, 
we can always choose a gauge transformation accordingly such that 

~^(0),ij = ~^(P),ij ~ ^ikl^k^lj — Texcess, ij ■ (66) 

We assume, therefore, that T( ) is symmetric and identical to T excess from now on. 

The gauge fixing can be performed because the following redefinition of T in the 
second integral of Equation (|58jl , 



a/ 



TijSRj = l {0)ij 5Rj - e m V k {VijSRj) , (67) 

does not change the value of the integral. Note that ^ = T(i)^ — eikiVijV k^Rj is 
also modified by the gauge transformation. In this context, we would like to point out 
also that a more general transformation of the tensor operator, 



J^R, =J ij 6R J +e ikl V k (w lj 6R^ , (68) 



where 



W 4J = Yl W(n),yfci-*n Vfc ■ ■ ■ V fc „ , (69) 

n=0 

is itself another tensor operator, also leaves the free energy unchanged. Equation (jfiHj) 
will be used later. 

Equation ()58|) will provide the connection between the variational quantities defined 
in Equation (JSTtj) and the mechanical quantities developed in the effective (Gibbs) 
description of the membrane mechanics, once a functional dependence of the cell 
deformation SR on the deformation of the dividing surface, SR, can be established. 
In the following two Subsections, we consider situations where the cell deformation 
SR is a local functional of the deformation of the dividing surface, SR, together with 
the shape itself, R, and possibly other fields like the density fields ha- We will take 
the parametrization of a material "particle" in the membrane-fluid system to be the 
(£\£ 2 ) /^-coordinates it had before the deformation took place. Given the requirement 
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that the mechanics of the system be invariant under rigid translations, we may propose 
a general form for a description of the cell deformation: 

SR = S(R + hn) + E ■ dJR 

= SR-ht a (n-d a SR) +E ■ dJR . (70) 

Each E for a = 1,2 is again a tensor operator. It should be clear to the reader that 
the differential operator d a is included in the above equation to ensure that SR = 5R 
under a rigid translation. 

The specific form of E depends on the particular model under consideration, and a 
physical example will be given at the end of this Section. A couple of general properties 
of E can, however, be readily derived. The first is a consequence of the fact that under 
a rigid, infinitesimal rotation of the system, SR = £ x R and SR = $, x R, where £ is 
a constant vector. Applying this requirement to Equation (f7T)|) leads to a condition on 

E* x t a = . (71) 
The second property, 

- = , (72) 

h=0 

is an obvious one, an expression of the fact that SR should be equal to SR at the 
dividing surface where h = 0. 

4-2. Translational shape deformations 

In this Subsection, we assume that SR represents a deformation resulted from a local 
translational movement of the dividing surface. Inserting Equation (|T0|) into the general 
functional form of the cell free energy, Equation (|58j) and performing the integral over 
h, we get the expression of the free energy in terms of surface-related quantities only, 
which reads as 

°Ft .excess = ' dA 



- f -SR-(nx t^) ■ t a (n ■ dJR) 

- t° ■ dJR + D a (T a -6R + f°- SR^j . (73) 
Two vector operators, T and are involved in the above Equation. T is given by 

f S / ffl (.-M + ^/„.r, (74) 
and is a tensor operator, but contains a regular, non-operator term. 

CO 



n=l 



is a full differential operator which is related to T. We are not, however, particularly 
interested in the operator , and do not, therefore, provide its explicit expression. 
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r • dR can always be reorganized into the following form, 

f a ■ d a 5R = 7? ) • 9JR + D a (7« ■ 8R) , (76) 
which defines a regular vector 7"^ and a vector operator 

oo 

n=l 

As we will shortly see, these two quantities will add terms to the stress tensors T a and 
Q a that has not been included in the previous formalisms P] QUI |2U], and they will 
be discussed further in a specific example later. Using Equation (|7fi|) we can rewrite 
Equation (J73J) as 

^S,exccss = J dA{-[f-D a (7f 0) + (n x r int ) ■ t Q n)] • 

(T a - 7 f 0) - (n x Tint) • t a n) ■ 6R 



+ D C 



+ [t:-^).5R 



(7f 



Comparing the above equation with Equation we can relate the excess 

mechanical quantities appearing here to the variational quantities / rs and S? \ derived 
from Equation (|56|) . Immediately we have 

T a = Sf 0) + 7f 0) + (n x r int ) ■ fn , (79) 

and 

f = f a + D a [ 7 f 0) + (n x r int ) ■ fn] . (80) 

The above two equations of connection constitute a new result, which has not 
been discussed in the existing literature so far. It is clear from the equations that the 
functional form of / alone is not sufficient for deriving the excess linear stress T a and the 
surface density of the resultant force, /: In addition we need the knowledge of 7^ and 
Tint, which are non-zero in general when external forces act on the membrane at positions 
that are not precisely on the dividing surface. Only in the case where there exist 
no such spatially distributed external force, both 7"^ and Tj nt vanish. Consequently, 
Equation (|79j) reduces to S?^ = T a , and we recover a connection used canonically 
earlier [TCI 120] ■ Similarly, Equation (JHUj) reduces to f = f rs . 

It may help to put the discussion in this Section in a context if we briefly mention an 
application of this variational approach: the establishment of equations of mechanical 
equilibrium for a membrane, or, the "equilibrium-shape" equations as they are called in 
the literature fZlEHHj- Once T a and / have been determined based on the derivation of 
Equation (|79*jl and Equation (jHUj) both from a given functional of the membrane excess 
free energy, /, and from given knowledge of 7^ and Ti nt , the equation of mechanical 
equilibrium that corresponds to force balance can be established from Equation (jSOj) to 
be 

/ + /■* + »• (T + -T-)| r=jR = 0, (81) 
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where 

/_ = /*(!- 2M, + *K) (/^ - - rL) , (82) 

— — J- 

with / rea i jext being the external force in the microscopic system and / cxt the external 
forces that are assigned to the bulk fluids in the Gibbs system. A simple example of such 
an application is given, for instance, by a membrane system, the mechanics of which is 
described entirely by the Helfrich bending-energy functional [Sj and where there is no 
external force applied. In this case, the relevant equilibrium shape equation is set up by 
setting 7"^ and T int equal to zero and by equating the component of / normal to the 
dividing surface to the "pressure difference" across the dividing surface ||. The equation 
resulted in agrees with that obtained from the variational procedure presented in, for 
example, Ref. [T7j . 

4-3. Rotational shape deformations 

The force / and the linear stress T a emerge naturally from Equation (|78jl. which 
results from the variation in the cell free energy with respect to a translational shape 
deformation, 8R, of the dividing surface. In this Section we demonstrate that in a similar 
fashion, the torque r and the angular stress fl a will emerge from the variation in the 
same free energy associated with a shape deformation that results from infinitesimal 
local rotations of the dividing surface. 

An arbitrary rotational deformation of the dividing surface is described by 8R = 
£ X R, where £ = £ 2 ) is a vector field prescribing the extent of the local rotations. 
Substituting this into Equation (J70j) leads to the following expression for the whole-cell 
deformation: 

8R = £x R-ht a (Rxn)- d a £ + E* • d a (£ x R) . (83) 

We may now replace 8R in Equation ()58|) with the above expression and derive an 

effective expression for the excess mechanical free energy contained in the cell in terms 

of the surface-related quantities only. 

Before doing that, however, we first make a point that will simplify the calculation 

of the free-energy variation based on Equation (|58J) . We have argued already that the 

mathematical expression of the free-energy variation may involve a non-zero Tm^-fc in 

the tensor operator T. But this term does not appear explicitly in the identification 

of the linear stress and the force. It turns out that Tm ^ can be neglected in the 

identification of the angular stress and the torque also. To make this point explicitly, 

let's consider an infinitesimal rigid rotation of the whole system, i.e., 8R = C x R, 

under which the free energy of the cell remains invariant in the absence of any external 

|| Note that the term "pressure difference" should be understood correctly as that defined in the Gibbs 
description. 
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force. Inserting this expression into Equation (|58|) and using the fact that the linear 
stress tensor T excess is symmetric, we find 



= ^Xl .excess = C t J dV V; {j (^^ijk^jlk) ■ 



J4) 



Since C and E are arbitrary the divergence in the integrand must vanish. The same 
mathematical theorem that led to Equation (|fi5|) allows us to conclude that 

^(l),ijk e jlk — (^ijk^kl) , (85) 

where V is a mathematically well-defined tensor. Recollecting Equation (|68j). which 
expresses the gauge freedom, or the arbitrariness, in the definition of the tensor operator, 
we have 

~^{l),ijk = T(l),yA; + e ikl^{0),lj + Qm/ V m W( 1 ) i /j fc . (86) 

We can conclude then that an appropriate \N(\)jjk can be chosen such that 

T{i),ijk e jik = . (87) 

We will, therefore, assume that Tm^^e,^ vanishes in what follows. 

Having made the above point, we can now substitute Equation (fKHJ) into 
Equation (p)H|) and derive the effective expression of ^-Fg^xcess- Using Equation (|H7|) . 
performing the /i-integral and arranging the result in a revealing form, we arrive at 

Excess = - j dA [r - D a ^ 0) - D a [(n x Tfart) • r (R x n)] } • £ 



/ (is z/ Q 



-(nxT int )-t Q (Rxn) (88) 



Several new quantities are introduced in the above expression. In particular, the vector 
quantity, ^" )' an< ^ ^ ne vec tor operator, u>", are defined by the following operation of 
the vector operator, f , defined earlier in Equation (J74j) : 

f a -d a (H xR)=u,? Q yd a £ + D a (u>°-£) , (89) 

with 

oo 

*: = E<r^---^- ( 9 °) 

n=l 

Obviously, W/jU and <I>° are related to 7^ and 7™ defined earlier in Equation (|74|1. and 
for u;^ the connection is given by 

^(o)^X7fo)-7:xJi. (91) 

Another vector operator, fi", has a general form given by 

00 

^ = E n w" A %-%. (92) 

n=l 

and is related to T". 
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The effective angular stress Q a appearing in Equation (|88|) must now be related 
to the corresponding variational quantities associated with the functional of the surface 
excess free energy /. In order to do that, we still need to calculate, for the rotational 
deformation of the dividing surface, 5R = £ x R, the corresponding variation in /. 
Inserting the deformation expression into Equation (JHH)) and reorganizing the result, we 
obtain 

^ (y/gf) = - (R x /„) • £ + D a (Q Q ■ . (93) 
The vector operator Q has the general form, 

oo 
n=0 

Its first term is given by the following particular expression: 

oo 
n=0 

related to the S?^'"^ n, s already derived from Equation (J5fi|) and Equation (foTj) . 
Comparing the line integral in Equation (|HHj) with the complete-derivative term in 
Equation (J9*3"|) we find 

Q a = Qf 0) + wf 0) + (n x r int ) • t a (R x n) . (96) 

Finally, we can establish, for instance, the connection between the effective bending 
moments, N a , and the variational quantity, Q(o), 

N a = + wf 0) - R x (Sf 0) + Tfoj) , (97) 

where Equation (fTU|) and Equation (|75jl have been used. This equation may be written 
into a more revealing form, where N a is divided into two contributions, 

N a = AT - 7? x R . (98) 

The first contribution, N® s , defined as 

2V« ^ Q^ 0) - x 5f 0) = - (5 a - Sf 0) ) x R , (99) 

can be obtained solely from the free-energy variation. The second contribution, —7" x R, 
requires additional knowledge of the cell deformation and any non-zero external forces. 
Again, Equation (}98j) makes it clear that a description of the membrane mechanics 
based on a given free-energy density alone is in general not sufficient for identifying the 
effective mechanical quantities of the membrane. 



4-4- Example 

To make the so-far rather general and formalistic discussion more concrete, we end 
this Section with an example. In the example, the surface density of the excess free 
energy is a given function of temperature T, some number density fields n^'s, the 
mean curvature H and the Gaussian curvature K defined for the dividing surface, i.e., 
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/ = f(T, H, K, {tia})- A specific example of such a function may be found in Ref. |2l?] . 
A simple purely geometric example is given by the Helfrich free energy [3] 

/= ^(2H-C ) 2 + kK + ct, (100) 

where k, Cq, R and a are phenomenological constants which are called bending rigidity, 
spontaneous curvature, Gaussian rigidity and tension respectively. 

Following the framework we have already established, we first calculate the variation 
of / = f(T, H, K, {ua})- Using the following set of identities from differential geometry, 

HVg) =\v99 a ^ 9af s, (ioi) 

5H = - l -K^5g aP + \g af3 5K aP , (102) 

5K = - Kg^6g a/3 + L^5K a p , (103) 
and carrying out the variation, we obtain 

^ (V9f) = l<r aP 59ap - \ aP 5K a p , (104) 



where 



o al3 



\*P = --^Lg^ - %L L <# . (106) 

2dH y dK y ' 

Further using 

&9ap — t a ■ dfjdR + tp ■ d a 5R , (107) 
8K aP = n ■ DpdJR , (108) 
we can rewrite Equation (|104|) in its final form: 

-U (V9f) = - D a [{^ - A "*/) tp + Dp\ al3 n] ■ 5R 
V 9 

+ D a {[ (a^ - X^K^) t p + DpX^n 

- X al3 ndp] ■ 5R} . (109) 

Comparing the above equation with Equation (JHH)) we can read off the various variational 
quantities one by one. S? \, which will contribute to the linear stress T Q , is given by 

S? 0) = (a a/3 - X^K^) tf, + DpX al3 n . (110) 

It should be pointed out here that a similar expression was obtained in [SO] by a 
sophisticated approach where geometric quantities were treated as auxiliary variables, 
which were constrained to match the appropriate expressions in terms of the shape field 
R through the Lagrange multiplier formalism. Inserting Eq. ()105j) and Eq. (jlOfij) into 
Eq. (jllOJl we arrive at 

^ -si**. ' m ' 

^■"^'{-km^-m^)- (112) 
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/ rs , which will contribute to the surface density of the deformation-related force, /, is 
given by its components in the normal and the tangential directions, 

f rs n= D a (Sf 0) ■ n) + Sfa ■ t?K af3 

V OK dn A ) dH K ' 

-\ D ^m- L °" D ^M' (113) 



/rs -t a — Dp • t a j - ■ nK a/3 



= -VnA^-. (114) 
dn A 

The simplicity of Equation ()114j) can be traced back to the reparametrization invariance 
of the free energy [HI] • 

To identify the variational contribution, N^ s , to the bending moments, we find first 
from Equation ()1U9|) that 

= - X af3 n , (115) 

Substituting these two results into Equation (j99|l immediately leads to 

N« = n x (X a %) = X ap e p % . (117) 

The formulas above can straightforwardly be applied to the purely geometric 
Helfrich free energy in Equation (jlOOjl . The results match the corresponding expressions 
derived through the Nother approach in [2*Uj . 

Note that N^ s ■ n = in this case. In general, however, the 2V£j obtained from 
the direct variation of a given free-energy density does not a priori have a zero normal 
component. On the other hand, the condition N a ■ n = must be satisfied, as we 
have already discussed in Section We will see in the next section that the normal 
component of iV° is completely arbitrary, and thus we can always choose a gauge in 
which it vanishes. The same is true for the 7" x R part, and therefore the condition 
N a ■ n = can always be satisfied (as it should be). 

We have already pointed out in Section 0] that a free energy alone is not sufficient 
in general for the complete determination of the effective mechanical quantities. Besides 
Tint, at least two additional quantities, 7^ and w" ^, are needed, which bear information 
on the details of the three-dimensional deformation of the membrane-fluid system. To 
illustrate how those quantities can be calculated, we consider an example, where the 
membrane changes thickness when its shape is deformed. We will later show how 
this example can be applied to membranes with a certain kind of active proteins 
incorporated. Mathematically we can model a change of thickness by the three- 
dimensional deformation 

5R = 5 (R + hn) + 5hn , (118) 
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where the change in the transverse direction will be written as an expansion to second 
order in h 

Sh = -2 (Ci + CihH) h-^- + 2( 3 h 2 5H . (119) 

Ci, C2 and C3 are three dimensionless phenomeno logical parameters. The term involving 
Ci + CzhH describes the effect of a mechanism where the membrane thickness shrinks 
(if (1 + (,2hH > 0) when the area of the membrane expands. The factor of 2 is chosen 
in order to get agreement with the definition of the Poisson ratio in elasticity theory 
[32] • The term characterized by parameter (3 models a type of structural changes in the 
system when the membrane is bend. An example of the effects which could be described 
by such a term is the flexoelectric effect [33] . 

Comparing the above equation with the equation of definition for the tensor 
operator Z a , (|7UJ1 . we find first that it is given by 

E Q = -2 (Ci + CnhH) hnt a - 2( 3 h 2 K a/3 nt f3 + ( 3 h 2 nnD a . (120) 

Going through the steps prescribed in Equation (|71jl and Equation (J7fij). we arrive at 
the following identification of 7^ and <^? y 

7(o) = — 2 (CiPcxccss + (2HQ CXCCSS ) t a —( 3 Q excess K al3 ti3—( 3 D a Q exccss n , (121) 



and 



where 



(0) 



R x Tfo) " C 3 Q™e a % , (122) 



Pexcess = j dh h (l - 2hH + Y? K) n ■ / cxcess , (123) 

Qexcess = J dh h 2 (l - 2hH + h 2 K) 71 • / cxcess . (124) 

We also obtain the "external" contribution to the bending moment 

- 7 ? x R = wf 0) - R x 7 f 0) = -( 3 Q excess e a % , (125) 
which is only tangential. Thus we can determine the bending-moment tensor to be 

= -\(m + 2 C3Qexcess) <f* ~ ^L<* . (126) 

We can give an explicit example of what the parameters Ci, C2 and (3 could be, by 
assuming that the membrane can be modelled as consisting of a material that mimics 
an incompressible fluid. A fluid particle at distance h from the membrane will then 
move to a distance h = h + 5h during a deformation in such a way that the fluid volume 
between the particle and the dividing surface is conserved. Mathematically we can write 
this criterion of volume conservation as (see Equation ()26jl ) 



= 8 ( J dhy/g(l- 2hH + h 2 K) 



^(h-h 2 H+h 3 K) ) . (127) 
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If we isolate 8h in the above expression and expand to second order in h we find 

Sh = - (h + h 2 H) 5 -^r + h 2 5H + O ((/i 3 ) . (128) 
Comparing with Equation ()119j) we see that for this specific example we have 

= C2 = C3 = \ ■ (129) 

An application of the considerations above can be found for a model of membrane 
activity proposed in [33] to explain experiments showing an increase in membrane shape 
fluctuations upon activation of certain proteins in the membrane. In this model the 
membrane proteins where activated by feeding them with the right energy source, which 
they then consumed to perform a specific task. While performing the task the proteins 
were assumed to change configuration or move nearby material in such a way that the 
proteins constantly pushed on their surroundings. Mathematically this was modelled 
as a distribution of forces added to the force-balance of the bulk fluids surrounding the 
membrane. In the language of this paper such a force distribution can be included by 
assigning it to the external force density / rea i ext . A way to write down the model would 
then be 



iWxtM = / dA ! dhF &ct {i\e,h)n5\r-{R+hn)) 
Jm J 



(130) 



where the specific expression for F act used in [H3] was based on each protein contributing 
a force-dipole: 

Fact = (Fa(p T - P l ) + 2HF& + p 1 )) [5 (h - w 1)-6(h + w 1 )} . (131) 

Here p^ and p^ represent area densities of the active proteins in the membrane, with j 
and I indicating the two possible orientations of an asymmetric transmembrane protein. 
F a and F' a are constants representing the strength of the active forces and their curvature 
dependence, and w^- are constant lengths giving the distances from the membrane 
where the forces act. If we switch form the coordinate r to (^ 1 ,^ 2 ,/i) and use that 
/real ext has to be balanced by an equal but opposite / real we can restate Equation (|13U|) 
as 

f^-^-T^PM- < 132 > 

We can define a Gibbs system corresponding to this microscopic system by stating that 
in the Gibbs system the bulk fluids are free of external active forces, i.e. / =0. Then 
/excess = /real ioi h ^ and we get the simple expressions 

= -/ d " F -- < 133 > 

Qexcess = - dh h? F act . (134) 



The expression given for 7"^ in Equation ()121j) . when supplemented by Equations ()129j) . 
(j!33J) and (|134jh is exactly the expression for the active contribution to the membrane 
stress T" ct derived in in a different manner. Thus the considerations on the effect 
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of external forces distributed in the transverse direction to the membrane developed in 
this paper provides an alternative point of view on the mechanics of this type of active 
membranes. We would like to point out that this kind of membrane activity is also an 
example of how the Gibbs description can be very useful. This was shown in jSH] where 
the Gibbs description derived in was used to calculate the fluctuation spectrum of 
an active quasi-spherical vesicle. 

5. Arbitrariness in the effective stresses 

As mentioned already in Section 0] in connection with the free-energy variation, 
Equation (|5fij) . there exists a certain degree of arbitrariness in the definition of the 
variational contribution S%\ to the linear stress T a . Related to that, there is also 
arbitrariness in the variational contribution Qf )> defined in Equation ()93|) . to the 
angular stress fl a . This issue has been mentioned in passing in some of the earlier 
literature [20], and it plays an important role in the recent article jSZj, where it 
is proposed to find solutions of membrane shape equations by equating an effective 
membrane stress tensor with a null stress that represents exactly the arbitrariness 
discussed here. In this Section, we discuss how to eliminate part of the arbitrariness the 
membrane stress tensor based on the understanding we have obtained of the effective 
mechanical quantities from the microscopic perspective. Moreover we will also present 
a geometry-based interpretation of the remaining part of the arbitrariness. 

It is not difficult to see that the free-energy variation formulated in Equation (JoT)j) 
is unchanged under a "gauge" transformation of S of the following form: 



n=0 

Equation ()135j) corresponds to a redefinition of the variational contribution to the 
linear and angular stresses 




(135) 



where W is an arbitrary vector-operator with the general form 



oo 




(136) 




S 




S^ + e a %W i0) , 
Q(o) - e aP dp (W x fl) 



(137) 




where 




oo 




(139) 



n=l 
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Note that there is more than sufficient freedom in the definition of W(*) that allows us 
to consider as an arbitrary vector-function %. 

The key to a physically meaningful elimination of at least some of the arbitrariness 
represented by the above equations is given by the bending moments. There is indeed 
a physical requirement on the bending moment, N a , namely, N a ■ n — 0. This in 
turn implies that it is natural to demand that the variational contribution, N^ s , should 
satisfy N® s ■ n = 0, and we will see that we have enough gauge freedom to ensure this 
requirement. In other words, we may use the arbitrariness in the stress definition to 
ensure that the physical requirement be fulfilled. Another condition is the condition 
that the tangential components of the variational bending moment M r a / = — iV° • Ve^ 
should be symmetric. This condition is a natural one when the bending moments may 
be described as derivatives of the free energy with respect to the symmetric curvature 
K a p. But, it is not a general one, based on an understanding of a microscopic origin. 

To see what that amounts to explicitly, we write out the "gauge" transformation 
for Nf s , which reads as 

N'« = N? B + e a % x W ( o) + e^dp W w . (140) 

It should be clear that the second term on the right-hand side of this equation makes 
it quite easy to choose the appropriate "gauge" that satisfies the physical requirement. 
It can be checked easily that both iV'° • n = and a symmetric M'^f can be achieved 
simultaneously by, for example, the following gauge 

W (0) = « . n )t fi - -t^n, (141) 

W w = . (142) 

It turns out that the two conditions mentioned above do not eliminate the 
arbitrariness, or "gauge freedom," completely. The remaining gauge freedom is described 
by a transformation 

W (0) = l -e^ (d a A ■t p )n- e a ? (d a A ■ n) t p , (143) 

W w = A , (144) 

where A is an arbitrary vector-function. It connects a sequence of iV" s 's and M^'s, 
which all satisfy the two conditions. In terms of iV" s and M%f, the transformation reads 
as 

N'l = N% + (e^5\ + ^V 7 ) {d p A ■ t 7 ) t s , (145) 
or equivalently, 

M' r f = M?f - - (d a A ■ t 13 + d p A ■ t a ) + <9 7 A • . (146) 

% Readers with knowledge of de Rham Cohomology (see for instance [HH|) should be able to convince 
themselves that at least locally (|135H , 1|137|) and l|138|) constitute the most general form of arbitrariness 

— Q' 

in S • SR, Sv ) an< ^ Q(o)j provided that their divergences are well defined quantities. 
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5.1. Arbitrariness from the Codazzi-Mainardi equations and theorema egregium 

No more obvious conditions that are physically meaningful can be imposed on 
the effective stresses derived variationally to eliminate the remaining arbitrariness, 
represented by the vector function A. In this Subsection, we develop a geometrical 
interpretation of A. This interpretation may help a user of the membrane mechanics 
make the most judicious choice, or gauge, in practice. 

The geometrical interpretation is derived from the Codazzi-Mainardi equations, 
given in Equation (|A.16|) . and theorema egregium, formulated as in Equation (|A.17|) . 
To describe qualitatively why the arbitrariness should be related to these equations, 
we recall first that a surface is uniquely defined, modulo rotations, translations and 
reflections, by the two tensors g a p and K a/3 [39]. The metric tensor g af3 describes the 
local extension/compression of the surface, the type of deformation from which the 
linear stress originates; the curvature tensor K a p describes the local bending of the 
surface, apparently the type of deformation that is associated with the angular stress 
or the bending moments. Were g a p and K a p independent of each other, we would be 
able to derive a a/3 and X al3 , generally and uniquely, from Equation (jlU4|) as functional 
derivatives of the free energy, and then determine the stresses and bending moments 
uniquely, as in Equation and Equation (jll7J) . g a p and K a/3 can not, however, be 

varied independently, due to the fact that they are tied together precisely through the 
Codazzi-Mainardi equations and theorema egregium. In other words, general bending 
and stretching of a surface can not be decoupled completely. Thus, the geometrical 
constraints should be related to the arbitrariness of the stress tensor. 

To show the connection more explicitly, we note first that the Riemann tensor TZ a ^s 
for two-dimensional geometrical manifolds can be written as 

7^«/?7<5 = -^£af3£-y6 ■ (147) 

The full non-trivial content of theorema egregium is thus captured by Equation (jA.23)) : 
H = 2K . (148) 

We then consider a "contribution" to the free-energy that reads as 

f = u n (K/2-K) . (149) 

This term is identically zero, and, therefore, makes no real contribution to the total 
free energy. The contributions to the linear and angular stresses derived directly from 
the variation of this zero free energy, using the definitions in section |Appendix A| and 
equations (|104|) . ()111|) and ()112|) . are, however, not zero, although their divergences 
must be zero. 

To calculate the variation, we first evaluate the variation of the Christoffel symbols, 
which are given by 

= \a' S (Dfs6g Sa + DJg sp - D 5 5g a p) . (150) 
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From this we can then derive, based on Equation (jA.18j) . 

^\ sp = D 5 (srj-D,(5Tl 5 ) , (151) 

and in turn 

5K = (-yS Q/3 + D a D P ~ g a ^D 7 D^j 6g a(3 . (152) 
Using the above result, we finally arrive at a formal expression of the free-energy 



variation: 



-^=S {y/gf) = \o-^8g a p - X af3 SK a/3 



where 



+ D^ n g ap - D a oo n g^)6g afS ] , (153) 
\ aP =oj n L aP . (155) 



Note that the boundary term inside the square brackets in the last line of Equation (jn 
is invariant under both rigid translations and rotations, and will not, therefore, 
contribute to the stress tensor or the bending moment. 

We can now use Equation (|1(J7|) and Equation (|1(J8|) to read off from Equation (|153|) 
the contributions to the stresses, which are 

M r <f = \ aP = uo n L aP , (156) 
5f 0) • n = D j3 \ l3a = L^d p u n , (157) 
Sf« 0) • t p = a a ? - X^K^ = D a D (3 u n - g aP D^D^uj n . (158) 

It is conceptually obvious that the above contributions necessarily imply a certain 
degree of arbitrariness in the definitions of the stresses, since the starting point is a zero 
free energy. That this is the case can also be seen, if we make a gauge transformation 
of the form given in Equation ()143j) and Equation ()144j) with A = u n n: the stresses 
become zero in the new gauge. 

The Codazzi-Mainardi equations are related to the tangential components of the 
vector function A. The proof is similar to what has been given above. We give only a 
brief presentation. The non-trivial content of the equations is captured by 

DpK p a = D a K% . (159) 

We consider, therefore, the following zero free energy 

/ = u a (DpK p a - D a K%) . (160) 
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Performing the variation in the same way as before, we get 



M r <f = - {D a w p + D p w a ) - D^g aP , 

S a ■ n = D p M Pa , 

S? 0) ■ t p = K a D p w~< - K^ s D 7 w 5 g aP 



(161) 
(162) 




(163) 



A useful mathematical fact for doing these calculations is that in two dimensions, the 
following identity holds for any tensor 



These stresses in fact represent the rest of the remaining gauge freedom, since they 
become zero under the gauge transformation defined by A = uj a t a . 

6. Conclusion 

Canonically, mechanics of a fluid membrane is effectively described as the mechanics 
of an infinitely thin surface. In other words, the effective description must reflect the 
physics that is associated with the finite, albeit microscopic, thickness of the membrane. 
Moreover, the effective description is formulated either in terms of the concepts of 
mechanical deformation and stresses, or in terms of mechanical free-energy functions. 
The connection between these two descriptions has been a subject of discussion in the 
existing literature QjJ EH ED] - A number of issues in this context have not, however, 
been addressed clearly and fully, or at all. In particular, they include the issue of what 
the microscopic origins are of the effective mechanical stresses, as well as the issue of 
what the connection is between the two descriptions when external forces act on the 
membrane system under consideration. 

It is the main purpose of this paper to address these issues, in a way that both casts 
the different existing works in a coherent framework and extends the current framework. 
To this end, we have approached both of the effective descriptions from a microscopic 
perspective, where the fluid membrane is treated as a microscopically thin layer, 
with highly inhomogeneous material and force distributions in its transverse direction, 
and have demonstrated unambiguously how the effective descriptions arise from the 
microscopic perspective, i.e., the microscopic origins of the surface mechanical stresses. 
Moreover and more importantly, we have, facilitated by the microscopic perspective, 
established a general connection between the mechanical-stress based and the free-energy 
based descriptions. Naturally, in doing so we have recovered the canonical connection 
already established in the existing literature, but we have also addressed a specific issue 
that is important to practical applications of the energy-based description of membrane 
mechanics: the issue of arbitrariness involved in identifying the mechanical stresses from 
a given free energy. We hope that we have provided more insight into the issue and 
that the understanding will facilitate the utility of both of the descriptions of membrane 



(164) 
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mechanics. Furthermore, we have worked out the connection for situations where there 
act external forces on the system. This is a non-trivial extension to the existing theories. 
We expect that it will be relevant and useful to further studies of membrane mechanics, 
both theoretically and experimentally. 
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Appendix A. Differential geometry of surfaces 

In this appendix we briefly review the mathematical language of differential geometry 
of two dimensional surfaces. A more comprehensive introduction can be found in Refs. 
[HUES], for example. 

The shape of a two-dimensional surface is represented by a space-vector function 
R = R^ 1 ,^ 2 ). The variables C, 1 and £ 2 are internal coordinates corresponding to a 
parametrization of the surface. At each point on the surface a basis for three dimensional 
vectors can be established. Two of them are tangential vectors, defined as 

dR 

t a = d a R= — , (A.l) 

where a = 1,2, and the third is a unit vector normal to the surface, given by 

n = . (A.2) 

|*i x t 2 \ v ' 

Local geometry of the surface is characterized by two surface tensors, the metric 

tensor and the curvature tensor. The local metric tensor is defined by 

g a p = t a -tp . (A. 3) 

It has an inverse, g a/3 , which satisfies, by definition, 

2 Q % 7 = ^ , (A.4) 

where 8" is the Kronecker delta and where the repeated Greek superscript-subscript 
indices imply summation following the Einstein summation convention. The metric 
tensor and its inverse are used to raise and lower Greek indices as in the following 
example: 

t a = g a % , t a = g a ^ . (A.5) 

The curvature tensor K a/3 is defined via the second derivatives of the surface shape 
function, 

K aj3 = n ■ d a dpR . (A. 6) 
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From it the scalar mean curvature H and the Gaussian curvature K can be obtained: 

H = \g aP K a0 , (A.7) 

K = det g ap K Pl . (A.8) 

Two other tensors will also be introduced here, 

e a p = e aP V9 , ^ = ^jy/g , (A.9) 

where e Q/ g = e a/3 with e n = e 22 = and e i2 = — e 2 i = 1 are tensor densities, and 
g = det g aj 3 is the determinant of the metric tensor. 

Expressions of covariant /contravariant differentiations of vector and tensor 
functions defined on the surfaces are facilitated by the use of the Christoffel symbols, T^. 
One instance, which will become particularly useful later, is the covariant differentiation 
of a surface vector function, w = w a t a , given by 

D a w p = d a w p + w^ ia . (A. 10) 

The Christoffel symbols can also be defined as certain combinations of the derivatives 
of the metric tensor, namely, 

r I/3 = (dagps + dpg 5a - d 5 gp a ) . (A. 11) 

It follows that the covariant divergence of w a can be written as 

D a w a = ^-d a (^w a ) . (A.12) 

The area of a local differential element of the surface is given by 

dA= y /gd£ 1 d£ 2 , (A. 13) 

an expression which will be repeatedly used in surface integrals. 

Two equations that are frequently used in this paper are the Gauss formula 

D a t p = K aP n , (A. 14) 

and the Weingarten equations 

D a n = -K af3 tf . (A. 15) 

Equation (|A.14|) follows from the definition of the curvature tensor and the definition of 
covariant differentiation, and Equation (|A.15|) is obtained by differentiating n ■ n = 1 
and n • t a = and solving for D a n. 

We will also from time to time weed both the Codazzi-Mainardi equations 

D a K Pl -D p K ai = Q, (A. 16) 

and an equation related to the famous theorema egregium by Gauss 

Tla^S — K al K/3S + KasK/3^ = , (A. 17) 

where the Riemannian curvature TV ' aS/3 can be calculated from the Christoffel symbols 
by 

We*? = d sKp - dpTi, + r] s r a , - r eP r aS . (A. is) 
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Another way to define the Riemannian curvature is through covariant differentiation 
of an arbitrary vector t> 7 : 

(D a D p - DpD a ) v* = K\ ap v 5 . (A. 19) 

In other words TV ' aS/3 measures the degree to which covariant differentiations commute. 
Contracting the Riemannian curvature tensor we get the Ricci curvature tensor 

*W = , (A.20) 

and further taking the trace we get the Ricci scalar 

K = g^n aP . (A.21) 
If we contract Equation ()A.17j) twice with the metric and use the following identity 

K ai K\ = 2HK aP - Kg aP , (A.22) 
we see that the Ricci scalar and the Gaussian curvature are intimately related as 

7l = 2K . (A.23) 

This is, in fact, one way to state theorema egregium. 

Theorema egregium and the Codazzi-Mainardi equations relate the intrinsic 
geometrical properties of the surface, i.e. those associated with the metric g a/ 3, to 
the extrinsic properties associated with K a p. In the main part of the paper it is shown 
how these relations lead to some arbitrariness in the definition of stresses in membranes. 

There is a symmetric tensor which will be used sufficiently often to merit a symbol 
on its own. It is defined as 

= e^e^Kse = 2Hg^ - . (A. 24) 

This tensor is proportional to the inverse of the extrinsic curvature tensor (when the 
inverse exists) 

L^Kjs = K5 P S . (A.25) 
It can be seen that 

DpL^ = , (A. 26) 

which follows from 

= , (A. 27) 

and from Equation ()A.16j) . 
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Figure Bl. A triangle bounded by coordinate curves and the curve c. 



Appendix B. The form of the membrane stress tensor 



In this appendix we will establish the form of the stress tensor given in Equation (J2J). 

To establish Equation (J2J) let us consider a small triangle in the membrane surface, 
as illustrated in Figure IB11 The triangle is bounded by a curve c and two other 
coordinate curves crossing each other at a point P. The outward pointing unit normal 
to the coordinate curve is u^ l > = t\ x nl Jg\\. The edge of the triangle formed by this 
curve has a length ds(i) = y/gTid!; 1 = \[g^ 22 v 2 ds, where d£} is the coordinate length of the 
curve and ds is the length of the part of curve c that contributes to the triangle. Similarly, 
the outward pointing unit normal to the other coordinate curve is = n x t 2 j ' y/g<n, 
and the length of the relevant part of the curve is ds( 2 ) = y/g-ndS, 2 = \[g^v\ds. If we 
define 



T 1 = -Vg^T {ul2)) ,T 2 = - v9®T (I/(1)) , (B.l) 

then the force on the £ x coordinate curve is (— T 2 / ^/g 22 )ds^, and the force on the £ 2 
curve is (— T 1 / g ll )ds{2)- The sum of the forces on the triangle is 



T(y)ds 



-a ds (])_ 

22 



,11 



(T (l/) - T 2 v 2 - TV) ds . 



(B.2) 



The sum of the forces should decrease at least proportionally to the area of the triangle, 
as the triangle is reduced by moving curve c closer and closer to Point P, in order that a 
force density of normal magnitude exist. In other words, [T(y\ — T 2 u 2 — T 1 ^) should 
approach zero as ds approaches zero, from which Equation (J2J follows. 
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